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ABSTRACT 

■ Aims. We present an analysis of the accreting X-ray pulsar system 4U 1907-1-09 based on INTEGRAL data. The main focus of this analysis is 
a study of the timing behavior of this source. In addition we also show an analysis of the 5-90 keV spectrum. 

Methods. The data were extracted using the official INTEGRAL software OSA 5.1. Timing analysis was performed using epoch folding and 
pulsar pulse phasing. 

Results. We have measured 12 individual pulse periods for the years 2003 to 2005. We confirm earlier RXTE results that during 2003 the 
. spin down became slower and show furthermore that after this phase 4U 1907-1-09 started to spin up with Ppuix = -0.158 s yr"' in 2004. The 
' similarity of the pulse period histories of 4U 1907-1-09 and 4U 1626-26 suggests that accretion onto an oblique rotator, as recently proposed 
, by Pema et al., is a possible explanation for this change. 
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1. Introduction 



The wind-a ccreting High Mas s X-ray Binary system 
4U 1907-1-09 JGiacconi et al.ll97ll) cons ists of a neutron star in 
an eccentric (e = 0.28) 8.3753 d orbit din 't Za nd et al."l998') 
around its companion, which has been identified optically 
with a highly reddened wy - 16.37 mag star Jschwartz et alJ 
[1980^. Using interstellar atomic lines of Na i and K L lcoxeTaLl 
(2005) set a lower limit of 5 kpc for the distance. According to 
this value a lower limit of the X-ray luminosity a b ove 1 keV 
is given by 210^'' erg s^' lin 't Zand et al.' 'l997l). 'Cox et ajJ 
(2005) al so confirm earlier suggestions (Schw artz et al. 1980t 
iMarshall & Rickettslll980t Ivan Kerkwiik et al.Hl989l) that the 
stellar companion is a 08-09 la supergiant with an eff'ective 
temperature of 30500 K, a radius of 26 Rq, a luminosity of 
5 10^ Lq, and a mass loss rate of 7 lO^^Moyr"'. Note that 
the presence of X-ray flaring seen twice per neutron star orbit 
( Marshall & Ricketts 1980, see also Fig. [0 had led some 
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auflio rs (e.g., iMakishima e7aDll984t Ilv3ll986l: ICook & Pag3 
11987") to the suggestion of a Be star companion, however, this 
classification would require a distance of < 1 .5 kpc, which is in 
contradiction to the sig nificant interstellar extinct ion measured 
in optical observations Jvan Kerkwiik et alJl989l) . 



Similar to other accreting neutron stars, the X-ray con- 
tinuum of 4U 1907 H-09 can be described by a power- 
law spectruni with an exponentia l turnover at ~13key 



dMihar. 




lin 'tZand et 


al.lll997l ICusumano et al.lll998h. 


The spectrum is modified by str 
tion with a column A^h - 1.5-5.7 


or 
1 


ig photoelectric absorp- 
(fi^ cmT^ ( Schwartz et alJ 


1980; 


Marshall & Ricketts 


1980 


Makishima et all Il984 


Cook & Page 1987;. Chitnis et al.J993;.Cusumano et al.ll998l). 



As the absorbing material is situated in the dense stellar wind, 
A^H is strongly variable over the orbit. The column is max- 
imal between the end of the primary X-ray flare and the 
start of the secondar y flare (Fig. ^ of the orbital light curve 
dRoberts et alJ 1200 ll) . Fluorescence in the absorbing mate- 
rial gives rise to an Fe Kg line at 6.4 keV with an equiva- 
lent width of ~60eV dCusumano et alJlTggsl) . At higher ener- 
gies, the spectrum exhibits cyclotron resonant scattering fea- 
tures (CRSF) at ~19keV and ~40ke V dMakishima et alJl992t 
lMiharall995tlCusumano et alJl998l). 
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With a pulse period of ~ 440 s, 4U 1907+09 is a slowly 
rotating neutron star. Si nce the discovery of the pulsations by 
iMakishima et alJ ( 1 19841) the neutron star has exhibited a steady 
linear spin down with an average of Ppuise - +0 .225 s yr~' from 
Pp„i.e = 437.5 s in 1983 to 440.76 s in 1998 (iiiT't Zand et all 
\l9m . Recently, iBavkal et al.1 tOQ^ reported a decrease in 
^'puise, which in 2002 was ~0.115syr ' and therefore ~0.5 
times the long term value. 

In this paper we report on an analysis of INTEGRAL ob- 
servations of 4U 1907+09 performed since 2003. In Sect.|2]we 
describe the INTEGRAL data reduction, followed by a discus- 
sion of the broad-band X-ray spectrum of the source (Sect.|3j- 
We then present the results of the long term X-ray timing anal- 
ysis (Sect. im i. showing that the source underwent a torque re- 
versal towards a spin up (Sect. R3l . We summarize and discuss 
our results in Sect. 13 

2. Observation and data reduction 

In this paper we present the analysis of INTEGRAL ob- 
servations of 4U 1907+09. The International Gamma-Ray 
Astrophysics Laboratory (INTEGRAL; Winkler etal. 2003) 
was launched in Octo ber 2002. The two main instruments 
IBIS (1 5 keV-lOM eV: 'Ubertini et allExB ) and SPI (20keV- 
8 MeV: fVedrenne e t al. 2003) are supplemented in the X-rays 
by an auxiliary instrument with a s maller field of view, the X- 
ray monitor JEM-X (3 keV-35 keVi lLund etaL>2003i) . For our 
analysis we used data from IBIS and JEM-X. 

The Imager on Board the INTEGRAL Satellite (IBIS) 
is a coded mask telescope with a fully coded field of view 
(FCFOV) of 9° X 9° and 1 2:9 full width at half maximum 
(FWH M) angular resolution JUbertini et al .'2003: Brand t et alJ 
l2003h . Here, we used data from its upper layer, the INTEGRAL 
Soft Gamma-Ray Imager (ISGRI), which covers an energy 
range fro m 15keV-l MeV wi t h an energy reso lution of ~8% 
at 60keV dLebrun et alJl2003t iGros et alJl2003i) . 4U 1907+09 
is detected below ~90keV. The Joint European X-ray Monitor 
(JEM-X) consists of two identical coded mask instruments 
with a FCFOV of 4^8 diameter a nd 3. '35 angular resolution 
(FWHM), each jLund et alJl2003l) . Due to erosion of the mi- 
crostrip anodes inside the JEM-X detectors which leads to a 
loss in sensitivity, only one of the two JEM-X detectors is op- 
erating at any given time while the other is switched off. 

For our analysis we took into account all public data up 
to INTEGRAL orbit (revolution) 250 as well as data from 
the INTEGRAL Galactic Plane Scans (GPS) and the Galactic 
Central Radian Deep Exposure (GCDE). In addition we also 
used private data from a monitoring campaign on the Galactic 
microquasar GRS 1915+105 (PI J. Rodriguez). To avoid sys- 
tematic effects that start to become important for larger off- 
axis angles, we constrained the maximum off-axis angle be- 
tween the satellite's optical axis and 4U 1907+09 to 4?5 for 
IBIS and 2°.4 for JEM-X according to the respective FCFOVs 
of the instruments. In total, the analyzed data are spread over 
almost three years with a total on source time of ~2280ksec 
for IBIS. The log of the observations is given in Table |3 We 
used the Off-line Scientific Analysis Software, 0SA5.1, in our 
analysis. In the analysis of the coded mask data we also took 
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Orbital-Phase (T,„ = 0) 

Fig. 1. Upper Panel: Rossi X-ray Timing Explorer All Sky 
Monitor (ASM) 2-12 keV light curve of 4U 1907+09 folded 
on the orbital period of 8.3753 days. Lower Panel: Same for 
the IBIS 20-40 keV light curve. 

into account the presence of GRS 1915+105, 4U 1901+03, 
4U 1909+07, and IGR J19140+0951, by forcing the OSA to 
consider 4U 1907+09 and these four sources. 

Background subtracted light curves were obtained using 
ii .light, a tool distributed with the OSA to generate high 
resolution light curves. While the standard OSA tools build 
lightcurves by deconvolving shadowgrams for each requested 
energy band and time bin to obtain the light curves, which is 
only possible if the signal to noise ratio in the shadowgram is 
high, ii_light uses the Pixel Illuminated Fraction (PIE) to 
create light curves. For 4U 1907+09, this approach allows us 
to extract light curves with a resolution of 1 s. The count rates 
found by ii_light show no systematic difference between the 
single pointing s. Tests based on ob servations of the Crab neb- 
ula and pulsar ( IKrevkenbohm et abi2006.) and the fact that the 
count rates of 4U 1907+09 are not affected by the very vari- 
able flux of the nearby source GRS 1915+105 (Fig.|5J confirm 
the stability of our extraction method. The light curves were 
barycentered and corrected for the orbital motion o f the neutron 
star using the ephemeris of lin 't Zand et alJ ( Il998l) . as given in 
Tabled Orbital phase is at Tgo. 

3. Spectral Analysis 

For the spectral analysis we applied systematic eiTors of 1% 
and 5% for IBIS and JEM-X2, respectively, and used data 
taken in the energy ranges 18-90keV and 5-20 keV. Because 
of the switch over from JEM-X2 to JEM-X 1 after revolution 
171, we only used data taken before that revolution to avoid 
cross-calibration uncertainties between the two JEM-X detec- 
tors. This resulted in a total exposure time of 200 ksec for JEM- 



Table 2. Log of observations. 
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Observation 


Revolution 


Date 


On Source Time (IBIS) 




Mean Count Rate 






MJD 


ksec 




cps in 20-40 keV 


1 


48 


52704.1-52705.4 


100.61 


0.78-0.93 


1.26±0.10 


2 


57 


52731.1-52732.3 


72.12 


0.00-0.14 


2.44±0.13 


3 


59 


52738.3-52739.6 


107.13 


0.86-0.02 


1.66±0.11 


4 


62 


52746.7-52747.9 


72.78 


0.86-0.01 


2.76+0.13 


5 


67 


52762.9-52763.4 


33.48 


0.80-0.86 


2.12+0.19 


6 


68 


52764.6-52766.5 


61.66 


0.00-0.22 


1.92+0.14 


7 


69 


52767.8-52769.5 


107.41 


0.38-0.59 


1.66+0.11 


8 


70 


52771.2-52771.7 


32.71 


0.79-0.85 


1.87+0.19 


9 


135 


52965.7-52966.9 


104.37 


0.01-0.16 


1.50+0.12 


10 


172 


53075.9-53076.2 


23.03 


0.17-0.21 


1.39+0.23 


11 


174 


53082.6-53083.0 


23.09 


0.97-0.01 


2.26+0.24 


12 


176 


53089.1-53089.6 


30.83 


0.74-0.81 


1.07+0.20 


13 


177 


53091.0-53092.0 


51.80 


0.97-0.09 


1.30+0.15 


14 


186 


53117.0-53117.7 


37.18 


0.07-0.16 


2.37+0.18 


15 


187 


53120.7-53122.0 


65.47 


0.52-0.67 


3.70+0.14 


16 


188 


53123.2-53125.4 


75.72 


0.82-0.08 


2.55+0.13 


17 


189 


53126.3-53128.5 


62.89 


0.19-0.45 


1.46+0.14 


18 


193 


53138.1-53140.5 


185.49 


0.60-0.88 


3.60+0.09 


19 


231 


53253.4-53254.0 


33.69 


0.36-0.43 


1.32+0.19 


20 


242 


53285.4-53285.6 


19.16 


0.18-0.21 


2.02+0.26 


21 


243 


53288.9-53289.2 


21.24 


0.60-0.64 


2.27+0.24 


22 


246 


53296.4-53297.6 


100.89 


0.49-0.64 


2.79+0.11 


23 


248 


53303.9-53305.0 


144.48 


0.27-0.52 


1.78+0.09 


24 


249 


53305.7-53308.0 


141.95 


0.61-0.87 


1.54+0.09 


25 


250 


53309.1-53311.0 


100.83 


0.02-0.23 


0.97+0.12 


26 


255 


53324.3-53325.5 


98.74 


0.82-0.97 


1.60+0.11 


27 


295 


53443.0-53444.1 


99.40 


0.99-0.14 


1.37+0.11 


28 


305 


53472.9-53474.1 


106.13 


0.57-0.72 


3.11+0.11 


29 


315 


53503.6-53504.8 


72.15 


0.24-0.37 


1.88+0.17 


30 


361 


53640.4-53641.6 


91.67 


0.57-0.71 


3.16+0.11 




Table 1. Ephemeris used for the binary correction 
din 'tZand et all Il998l) . forb is the orbital period, Tgo the 
time of mean longitude 90°. All other symbols have their usual 
meanings. 



Parameter 


Value 




MJD 50134.76 




8.3753 d 


a sin i 


83 It-s 


Eccentricity 


0.28 


0) 


330° 



X2. For ISGRI we used all public data up to revolution 250 
resulting in an exposure time of 1258 ksec. We modeled the 
spectrum using an absorbed power-law, which is modified by 
the Fermi-Dirac cutoff (■Tanaka^l986) : 



exp 



-^cut E 
Efoid 



1 



(1) 



This continuum model alone could not describe the data ad- 
equately (x^^d - 3.10, see Tab. as strong absorption line 
like features at ~19keV and ~40keV remain (Fig.|3l. These 
features have been identified as cyclotron resonant scattering 
features (CRSFs) by .Makishima et al. (1992) and have also 



Fig. 2. ISGRI image of the 4U1907-H09 region in the 20- 
40keV band (exposure time about 100 ksec, data from revo- 
lution 246). 



been reported bv lMiharal (Il995h and ICusumano et al] ( Il998l) . 

The addition of a line at 36.0keV with a width of 3.7keV 
(Table 13 with a Gaussian optical depth profile (for definition. 
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see ICoburn et al.ll2002l) to the baseline continuum model re- 
sults in a improved fit (x^^^ - 2.49 for 47 degrees of free- 
dom). The F-test probability that this improvement is achieved 
just b y chance is 3.8 x 10"^ (see, however, Protassov et al. 
I2OO2I) . Modeling both features with Gaussians at 39.8 keV and 
18.9keV further improves significantly to x^si - l-OO- 
The F-Test probability for this case as opposed to the case 
without lines is 9.8 x 10"'", confirming the unambiguous de- 
tection of the cyclotron lines. For this model we also added 
the Fe-line wh ich has been seen, e.g., in the BeppoSAX data 
iCusumano etalJ ll998). The best fit value for the line energy 
is 7.1 keV. The width of the line could not be constrained by 
our data, so we fixed it to e V, corresponding to an equivalent 
width of 135 eV. 

To exclude the possibility that the CRSFs are due to incor- 
rect modeling of the continuum, we also used other continuum 
models' including the NPEX model, which consists of nega- 
tive and positive power laws with a common exponential cutoff 
factor fchara 1995): 



C(E) = Al +A2 ■ E^"^) exp 




Using an absorbed NPEX model, we obtain x^f.^ - 4.54 and 
also observe line like features at 19 and 40keV. Again we first 
fitted only one feature at 44.4keV, resulting in ^ X^^^i - 2.56. 
Fitting these features with two Gaussians at 42.0 keV and 
18.7keV improves the fit significantly resulting ^t^X^^A ~ ^-^^ 
(F-Test: 8.7 x lO"'^). See Table El for the complete set of 
model parameters. With the NPEX model it was not possi- 
ble to constrain the Fe-line, however the upper limit for a line 
with the same energy and sigma as in the Fermi-Dirac cutoff 
model is 9.7 10"'* photons cm"^ s"', consistent with the results 
obtained for the continuum with a Fermi-Dirac cutoff. The con- 
tinuum and line parameters measure d with INTEGRAL ar e 
consistent with the BeppoSAX results dCusumano et al.ll998l) . 
Our measurements also confirm the BeppoSAX result that in 
this source the second line is deeper than the fundamental 
line, a fact that can be explained by "photon spawning", i.e., 
Compton scattering with multiple photon emission where a 
photo n with energy 2Ecwc..\ results in two photons with energy 
iicyci dAlexander & Meszaroslll989l) . From our spectral fitting 
we obtain an average flux of 2.6 10-'"ergcm-^ s"' in the 20- 
40keV band. 

4. Timing Analysis 

4.1. The X-ray light curve of 4L1 1907+09 

The X-ray light curve of 4U 1907+09 shows a clearly pulsed 
signal with a period of ~441 s. The mean count rates in the 20- 
40keV band for our observations are shown in Table |2 The 
variations in count rate largely reflect the orbital phase depen- 
dent variability (Fig.Q). 

4U 1907+09 is also known to exhib i t flares on a 
timescale of hours ( Maki shima et alJ Il984t lin 't Zand et alJ 

' We did not use the highecut model due to the artificial fea- 
tures introduced in the fit residuals by the abrupt onset of the cutoff 
jKretschmar et alll997tlkrevkenbohm et alll999h . 




5 10 50 

Energy [keV] 



Fig. 3. Fit to the data using the Fermi-Dirac cutoff model in- 
cluding 2 cyclotron lines. The bottom panels show the residu- 
als, from top to bottom the case without lines, with 1 line, and 
with 2 cyclotron lines. 

Table 3. Best fit parameters for the Fermi-Dirac cutoff model. 
Shown are the best-fit absorbing column, A^h, the photon in- 
dex, r, the flux at 1 keV, A, the cut-off energy, E^ut, the folding 
energy, Em^, and the cyclotron line parameters (energy, E^yc, 
depth Dcyc, and width, cTcyc), and the reduced and the num- 
ber of degrees of freedom, dof. All uncertainties are at the 90% 
level for one interesting parameter. 







fdcut 






without line 


with 1 line 


with 2 lines 


A'h [10^- cm-2] 




10!? 


0+' 

"-0 


r 


2 01+°" 


9 1+0.15 


" -0.09 


A 


21+0 13 


"•■'^-0.13 


o.iit« 


[keV] 


31+1 


37!^ 


301 


EfM [keV] 


5+' 




'-2 


£cyc,i [keV] 






18 9+°* 


-^cyc,l 






35+""' 

"■■'-^-0.09 


o"cyc.i [keV] 






3 0+' " 


Ecyci [keV] 






39.8!*;° 


Ocyc.2 




' -0.2 


0.9!°:^ 


a-cyc,2 [keV] 




-'• '-2.2 




fipe [keV] 






7 1+" 


o-Fe [eV] 






(fixed) 


Ape [10-*] 






6!? 


xVdof 


154.91/50 


116.82/47 


41.90/42 


^red 


3.10 


2.49 


1.00 



S. Fritz et al.: Torque reversal of 4U 1907+09 



5 



Table 4. Best fit parameters for the npex model, i.e., the photon 
index of the lower energy power law, ai (the positive power law 
index, a2, was fixed at 2.0), the power law normalizations Ai 
and A2 as defined in Eq. |2 and the folding energy Zifoid. All 
other symbols have the same meaning as in Table|3] 







npex 






without line 


with 1 line 


with z hnes 


Nh [10-2 cm--] 


1- 


41 




a, 




0.98t?2 




A, [10-2] 


r] 


9+7 




A2 [10-"] 




2+' 


61 


£foid [keV] 




6 21+"-"*' 


5.16!°^^ 


£cyc,i [keV] 






18.7!«| 


^cyc,l 






0.32!«;2o 


o"cyc.i [keV] 






2-8!ii 


£cyc.2 [keV] 




AA 4.+2-5 


42 0+^ ' 






1 07+0.37 
-0.29 


77+'' ''* 


cr,.yc,2 [keV] 




'•"-1.2 


7.8!^:^ 


ArVdof 


227.05/50 


120.46/47 


50.88/44 




4.54 


2.56 


1.16 



Table 5. Flares in the X-ray Ught curve of 4U 1907+09. The 
peak rates are quoted for a 441 s binned light curve in the 20- 
40keV band. 



IVIJD 

53138.4 53139.0 53139.6 53140.2 

F"T"T"T"T"T"T"T"T"T"T"T'n"n"n"n"n"n"n"n"'T"T"T"T"T"'l 

: MJD 




6 12 18 6 12 18 6 12 
2004 May 13 2004 May 14 2004 May 15 

Fig. 4. IBIS (ISGRI) 20^0 keV Hght curve of revolution 193. 
The bin time is 200 s. The inset shows a close up of the flare, 
the bin time in this case is 40 s. X-axis numbers are hours. 

MJD 

47000 50000 53000 

4421 ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 



Revolution Center of Flare Duration Peak Rates 

MJD sec IBISctss-' 

187 53121.65 7000 063 181)9 

188 53125.03 5000 0.03 12.11 
193 53138.21 7000 0.61 22.56 
305 53472.93 2500 0.57 13.88 



ll998HMnkeriee etal.l200lh . We observed four such flares with 
count rate increases of at least 5cr over its normal level (Fig.0}. 
Three of them, in revolutions 187, 193, and 305, are associated 
with the main peak in the orbital lightcurve (Fig.n* while the 
flare observed in revolution 188 is linked to the secondary peak. 
Table|5lshows an overview of the properties of these flar es. 

Du ring flares, lin 't Zand et al. (.1998) and Mukeriee et alJ 
( I2OOII) also reported the detection of transient 18.2 s and 14.4 s 
quasi-periodic oscillations (QPOs). We therefore considered to 
calculate dynamical power spectra for the flare light curves to 
search for similar oscillations, however, simulated light curves 
showed that our INTEGRAL data are not sensitive enough to 
detect a QPO at the reported levels. 

4.2. Pulsar Period and Pulse Profiles 

It is possible to track the evolution of the pulse period of 
4U 1907+09 for more than 20 years. Since the first mea- 
surements, the source showed a steady spi n down rate of 
^'puise = +0.225 syr-' din 'tZand et alJll998l) . A recent anal- 



436 1 , , \ , , , , \ , , , , \ , , , , \ , , , , \ ^ 

1985 1990 1995 2000 2005 

Fig. 5. Evolution of the 4U 1907+09 pulse period over the last 
20 years. Diamonds indicate the results obtained in this work, 
other s ymbols indicat e data from|Makishima et al. ( 1984, plus- 
signs l ICo ok & Page (1987, asterixes), Mihara (1995, trian- 
gles), E^rZande^ala(il99^ squares), and B avkal et aL(i,200lL 
2005, crosses). From MJD 5 1980 onwards 4U 1907+09 shows 
a clear deviation from the long term spin down trend with a 
complete reversal to a spin up after MJD 53131. The solid line 
shows the fit of four distinct episodes of different, but constant, 
P. 

ysis of RXTE data, however, taken during MJD 51980-52340 
showed a sp in down rate which is 0.5 times lower than the pre- 
vious value jBavkal et al1l2005h . Using our INTEGRAL data, 
we can extend the pulse period evolution to the years 2003- 
2005 (MJD 52739-53504). 
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To measure the pul se period of 4U 19 07+09 we first per- 
formed epoch folding llLeahv etalJll983h on each revolution 
separately. Using these periods we derived pulse profiles for 
each revolution using the respective period. The pulse profiles 
were then fitted by a sum of sine functions to determine the po- 
sition of the trailing edge (see for example Fig at phase 1.0) 
of the pulse to obtain absolute pulse arrival times and the differ- 
ences between the arrival times of adjacent revolutions. From 
these time differences the approximate number of pulse cycles 
between those absolute times was determined. For sufficiently 
well sampled observations, those values are close to an inte- 
ger number These integer numbers are then taken to count the 
cycles throughout those well sampled subsets, and a model can 
be constructed to calculate expected pulse arrival times (trailing 
edges). The observed times can then be compared with the cal- 
culated times and a polynomial fit to the observed times yields 
values for the period and its derivative. Care was taken that 
within subsets no miscounting was possible. Between subsets 
the time gaps are too large such that they cannot be bridged 
uniquely and no statement about the period development dur- 
ing the gaps can be made. Using this procedure we were able 
to derive 12 periods with sufficient accuracy (see Table|6j. 

Fig- 13 shows the long-time history of the period evolution 
based on all available data. The line indicates the historic spin 
down trend with Ppuise = +0.225 syr ' from in 't Zand et al. 



ill997l) . The recent periods obtained by iBavkal et al J (I2OO5I) 

show a clear deviation from the historic spin down trend. The 
periods obtained from INTEGRAL (shown as diamonds), not 
only confirm this change of the trend but show a complete trend 
reversal from the historic long term spin down trend to a spin 
up trend from MJD 53131 onwards (see Figs. |5] and |6}. 

We tried several models to fit this overall period evolution^. 
The description of the data by three distinct episodes with a 
linear spin down from MJD 45576 to MJD 50290 (Ppuise = 
+0.235+0.001 s yr"' ), a parabolic turnover from MJD 50290 to 
MJD 53022 of the form P(r) = Po+Ppuis,fl(T-To)+PpuU,.o(T- 
Tof/2 with To = MJD51321, Ppuke.o = +0.121 + 0.001 syr"! 
and Ppuise.o = -(7.743 + 0.002) lO"-' s yr"^ , and a linear spin up 
with Ppuise = -0.147 + 0.006 syr-' from MJD 53022 onwards 



Table 6. Period 
measurements ar 



measurements of 4U 1907+09. 
! taken from IMakishima et al. 



resulted in;i^^j,j = 8.74 for 34 degrees of freedom. Modeling the 
period evolution of 4U 1907+09 with four distinct episodes of 
different, but constant, P results in a slightly better O^^d ~ 
7.2 for 34 dof, see Figs.|5]and|6jl. The four episodes are the his- 
toric spin down with Ppuise = +0.230+0.001 syr^' (MJD 45576 
to MJD 50610), the turnover with Ppuke = +0.114+0.001 syr"' 
(MJD 50610 to MJD 52643) and Ppui.e = +0.026 + 0.003 s yr"' 
(MJD52643 to MJD53131), and finally the spin up with 
/'pulse = -0.158 + 0.007 syr-1. 

Using the derived periods, pulse profiles for all revolutions 
were determined (Fig.0 by folding the light curves with the 
best respective periods. No change in the shape of the pulse 
profile is detectable, and the INTEGRAL pulse profiles before 
and after the reversal of P are consistent with each other In 
addition, we also obtained energy resolved pulse profiles (see 




in 't Zand et al . (19981 198), Mukeriee et al 
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1995, 


. I 


12001, 



H istoric 
(11984 

M95), 
MOl), 



Date 


Instrument 


Pulse Period 


Reference 


MJD 




s 




45576.5 


Tenma 


437.438 ± 0.004 


M84 


45850.7 


EXOSAT 


437.649 ± 0.019 


C87 


48156.6 


Ginga 


439.19 + 0.02 


M95 


50134.8 


RXTE 


440.341 ± 0.014 


198 


50302.0 


IXAE 


440.53 ±0.01 


MOl 


50424.3 


RXTE 


440.4854 ±0.0109 


B05 


50440.4 


RXTE 


440.4877 ± 0.0085 


BOl 


50460.9 


RXTE 


440.5116 ± 0.0075 


B05 


50502.1 


RXTE 


440.5518 ±0.0053 


B05 


50547.1 


RXTE 


440.5681 ±0.0064 


B05 


50581.1 


RXTE 


440.5794 ± 0.0097 


B05 


50606.0 


RXTE 


440.6003 ±0.01 15 


B05 


50631.9 


RXTE 


440.6189 ±0.0089 


B05 


50665 


IXAE 


440.95 ± 0.01 


MOl 


50665.5 


RXTE 


440.6323 ± 0.0069 


B05 


50699.4 


RXTE 


440.6460 ± 0.0087 


B05 


50726.8 


RXTE 


440.6595 ±0.0105 


B05 


50754.1 


RXTE 


440.6785 ± 0.0088 


B05 


50782.5 


RXTE 


440.6910 ± 0.0097 


B05 


51021.9 


RXTE 


440.7045 ± 0.0032 


BOl 


51080.9 


RXTE 


440.7598 ± 0.0010 


BOl 


51993.8 


RXTE 


441.0484 ± 0.0072 


B05 


52016.8 


RXTE 


441.0583 ± 0.0071 


B05 


52061.5 


RXTE 


441.0595 ±0.0063 


B05 


52088.0 


RXTE 


441.0650 ±0.0063 


B05 


52117.4 


RXTE 


441.0821 ±0.0062 


B05 


52141.2 


RXTE 


441.0853 ±0.0082 


B05 


52191.4 


RXTE 


441.1067 ±0.0046 


B05 


52217.2 


RXTE 


441.1072 ± 0.0077 


B05 


52254.3 


RXTE 


441.1259 ± 0.0074 


B05 


52292.0 


RXTE 


441.1468 ±0.0065 


B05 


52328.8 


RXTE 


441.1353 ±0.0090 


B05 


52739.3 


INTEGRAL 


441.253 ±0.005 


this work 


52767.1 


INTEGRAL 


441.253 ±0.005 


this work 


53083.9 


INTEGRAL 


441.283 ±0.005 


this work 


53121.1 


INTEGRAL 


441.274 ±0.005 


this work 


53133.4 


INTEGRAL 


441.297 ±0.005 


this work 


53253.6 


INTEGRAL 


441.224 ±0.010 


this work 


53291.3 


INTEGRAL 


441.201 ±0.005 


this work 


53314.0 


INTEGRAL 


441.188 ±0.005 


this work 


53324.7 


INTEGRAL 


441.183 ±0.005 


this work 


53443.4 


INTEGRAL 


441.154 ±0.005 


this work 


53473.3 


INTEGRAL 


441.139 ±0.005 


this work 


53503.8 


INTEGRAL 


441.124 ±0.005 


this work 



^ Fo r this analysis we exc luded the two period measurements by 
IXAE JMukeriee et al.ll200lh as outliers. Including them resulted in 
X^^^ values increased by a factor of 5. 



Fig lUfor an example). Again, the evolution of the pulse profile 
from below 30 ke V to high energies is identical for all revolu- 
tions: the profile exhibits a single peak at all energies with only 
a small secondary hump, indicating the potential presence of a 
second peak at energies <20 keV, where a clear double peaked 
pulse profile was observed previously fcMakishima et al...l984t 
ICook & Pag3ll987t lin 't Zand et alJll998l) . At energies above 



S. Fritz et al.: Torque reversal of 4U 1907+09 



7 



52800 



MJD 

53100 



53400 



441.4 r 



441.3 - 



441.2 



441.1 



441.0 




MAMJ JASONDJFMAMJ JASONDJFMAMJ 
2003 2004 2005 

Fig. 6. Close up of the periods found in this work. Between 
MJD 52643 to MJD 53131 the derived pulse periods show a 
spin down of the source with Ppuke - +0.026 syr"'. After 
MJD 53131 the accretion torque changes sign and 4U 1907+09 
spins up with Ppuise = -0.158 syr The line indicates the lin- 
ear fit to the periods as described in the text. 



60 keV, the source gets too faint and no pulsed emission is ob- 
servable. 



5. Summary and Discussion 

Since 2003, 4U 1907+09 was observed frequently with 
INTEGRAL. In this paper, we presented the results of the tim- 
ing analysis of ~2280ksec of IBIS data and of the spectral anal- 
ysis of 1258 ksec of IBIS and 200ksec of JEM-X2 data. Two 
cyclotron lines at ~19keV and ~40keV are de tected, consis- 
tent with earli er Ginsa and BeppoSAX results ('Mihara' U995t 
ICusumano et al...l9981 Assuming a gravitational redshift of 
z = (1 - (2GM/rc^))-^'^ - 1 =0.31 (assuming M = I.4M0 
and r = 10^ cm), from these CRSF observations a B-field 
strength of 2.15 x 10'^ G can be infeiTed. Four flares in the 
X-ray lightcurves with durations between ~2500 s and ~7000 s 
were observed, three of which are associated with the primary 
flare in the orbital lightcurve of 4U 1907+09. 

In addition to these confirmations of earlier results, our 
INTEGRAL data show, for the first time, a clear spin up phase 
in 4U 1907+09 (Table|6land in Figs.|5land|6}. After almost 20 
years of constant spin down at a rate of Ppuise = +0.225 syr"'. 



Os- 



m 



with the available data being consistent with P 
2002 and 2003 first indication s for a decrease in th e magnitude 
of spin down were detected llBavkaletalJl2005h . As shown 
in Sect. 14.21 the new INTEGRAL data show a torque rever- 
sal from ~MJD 53131 onwards, with the source now exhibit- 
ing a spin up with a rate of P^uUe - -0.158 syr'. So far, the 
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Fig. 7. Pulse profiles for diff'erent revolutions in the 20^0 keV 
band. Panels a-e show pulse profiles obtained during the period 
of spin down, and panels f-o show spin up pulse profiles. 
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Fig. 8. Energy resolved pulse profiles for revolution 193 
(^puise = 441.224 s). 



INTEGRAL results appear to be consistent with P = s ' (the 
Icr upper limit for P during the spin up is -4 10"^ s"' yr"^). 

The traditional interpretation of spin up in X-ray binaries 
with a strongly magnetized, disk-accreting neutron star is that 
the accretion disk is truncated at the inner edge o f the disk by 
the rnagnetic field of the netitr on star (see, e.g., lOhosh et al.l 
Il977t iGhosh & La'mblll979albl and references therein). For 
prograde disks, the magnetic coupling between the accretion 
disk and the neutron star and the transfer of angular momen- 
tum from the accreted matter onto the neutron star will lead to 
a torque onto the neutron star resulting in spin up. The general 
expectation of these models is that the accretion disk is close 
to equilibrium, i.e., the Alfven or magnetospheric radius, r,,,, 
where the accretion flow couples to the magnetic field, is as- 
sumed to be close to the corotation radius, r^o, where the disk's 
Kepler frequency equals the rotational frequency of the neutron 
star. For accretion to occur, r,,, < r^o, as otherwise the onset of 
the centrifugal barrier (the "propeller eff ect") prevents the ma- 
terial from falling onto the neutron star jlllarionov & SunvaevI 

Most accreting neutron star systems, such as Her X-1 or 
Vela X-1, show a secular spin up or spin down trend onto 



which short, large magnitude episodes of spin down or spin up 
are superimposed jBildsten et alJl997l) . These shorter episodes 
are often interpreted to be due to short term variations in 
M, which give rise to short-term torque fluctuations. In con- 
trast to these systems, GX 1+4, 4U 1626-67, and 4U 1907+09 
showed decade-long phases where P did not change its sign. 
GX 1+4 was seen to spin up with a very short characteris- 
tic timescale, P/P ~ 40 years, from the early 1970s until the 
source flux dropped below detectability in the EXOSAT era 
of the early 1980s. From 1985 onwards the source re emer ged 
with a strong spin down jMakishima et alJll988t iMonv et alJ 
'1991^, 'Chakrabartv et al.lll997l ). Recent INTEGRAL observa- 
tions show that this spin down has continued until at least 2004 
jFerrigno et alJl200a) . A second source with a long phase of 
unchanged sign of P is the ultracompact low mass X-ray bi- 
nary 4U 1626-67, which was spinning up between 1977 and 
1989, followed by a torq ue reversal and an e xtended spin down 
phase starting in 1991 ( iChakrabartv et alj[r997fe . In contrast 
to GX 1+4, no luminosity drop was seen during the torque 
reversal and with \P/P\ ~ 5000 years the characteristic spin 
up and spin down timescales are significantly larger. Provided 
that the torque reversal of 4U 1907+09 described in this paper 
does not change back in the next years, in its overall period 
behavior, 4U 1907+09 resembles 4U 1626-67: both sources 
have \P/P\ ~ several thousands of years and no change in 
luminosity was observed during the torque reversal. On the 
other hand, however, 4U 1907+09 is a high mass X-ray bi- 
nary with a much longer orbital period, and thus shares sim- 
ilarities with the low mass X-ray binary GX 1 +4, and further- 
more 4U 1907+09 switched from spin down to spin up, while 
GX 1+4 and 4U 1626-67 switched from spin up to spin down. 
Common for the three sources is that the magnitude of \P/P\ is 
similar for the spin up and the spin down. 

With the simple magnetic torquing model outlined above, 
the long distinct episodes of roughly constant P and the magni- 
tude and sign of P are difficult to explain: As discussed above, 
for systems close to equilibrium it is expected that ~ r^o- In 
reality, however, this is no t the case. In the case of 4U 1907+09, 
it was first pointed out bv lin 't Zand et alJ ( Il998l) that the mag- 
netospheric radius inferred from the cyclotron line measure- 
ments is rm ~ 2400 km, while r^o ~ 12000 km. Moving 
out to Too would require a magnetic field of ~10''*G, which 
is two orders of magnitude larger than the magnetic field 
deduced from the observed cyclotron lines. Alternatively, M 
could be significantly lower than the value inferred from the 
X-ray luminosity of the system. Assuming the efficiency of 
accretion is comparable to other systems, however, a low M 
would imply a significantly smaller distance to 4U 1907+09 
than a llowed by its optical re ddening, which seems equally un- 
likelv. lin"'t Zand et aLN 19981) argued that the presence of quasi- 
periodic oscillations (QPOs) and the long term systematic trend 
in the spin down are strong indicators for the presence of a 
accretion disk in the sy stem with a small inner disk radius. 
In analogy to GX 1+4 JChakrabartv et al.llT997i and refer- 
ences therein), these authors show that a transient retrograde 
disk with a duty cycle of 1-5% could provide a sufficiently 
strong torque to explain the observed spin down. This torquing 
would coincide with a short term increase of source luminos- 
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ity. lin 't Zand et alJ ( Il998l) speculate that the X-ray flares in- 
dicate the presence of the retrograde disk, since a switch to 
a retrograde disk would imply an increase in X-ray luminos- 
ity. They point out, however, that the ~ 1000 s duration of the 
flares seen in their RXTE data is too small with respect to 
the required duty cycle for a retrograde disk. We note that the 
INTEGRAL results show flare durations of several ksec du- 
ration (Tab. |5}, which is more in line with the expected duty 
cycle. Furthermore, the folded lightcurve of 4U 1907-1-09 in- 
dicates a brightening of the source around phase 0.6 (Fig. 0. 
With O.lPorb duration, this flare is long enough to easily ac- 
commodate the required duty cycle of the retrograde disk and 
to explain the observed P. We note, however, that not all of the 
brightening seen in Fig.Qlis due to the flares identified by us as 
removing these flares from the data still results in a clear peak 
at phase 0.6 in the folded lightcurve. Part of the brightening 
could therefore be also due to an enhanced M during this phase 
of the elliptical orbit of the neutron star, although we cannot 
firmly exclude the presence of smaller flares which would not 
be picked up due to our conservative definition of a flare. 

While the model of a retrograde disk can explain the mag- 
nitude of P/P, it is more difficult to reconcile the torque re- 
versal found with INTEGRAL with this model. First of all, 
as discussed above, the four intervals of different P represent 
distinct episodes of different torques on the neutron star In 
the retrograde disk model, it is difficult to understand why 
the characteristic spin up and spin down timescales appear 
so similar. Furthermore, these three phases would then corre- 
spond to different duty cycles of the retrograde disk, resulting 
in a change in observed flux. Such a behavior has been ob- 
served, e.g., in GX 1+4, where torqu e and luminosity are cor - 
related during the spin down phase ( Chakrabartv et alJl 19971) . 
No such correlations are observable in 4U 1907 H-09, where the 
flux has not appreciably changed during t he INTEGRAL an d 
the pointed RXTE observations (see also lBavkal et al.l l2005). 
In addition, an analysis of the orbit-averaged 2-12keV RXTE 
ASM lightcurves also does not reveal a change in the soft X- 
ray behavior of the source, and neither is a change in the X-ray 
light curve with orbital phase observed. 

The INTEGRAL observations also rule out the applicat ion 
of the torque reversals put forward bv iMurrav et alJ 11999) to 
explain GX 1 +4. In this model, the accretion disk is assumed 
to consist of several rings with opposite rotation. Between 
these rings a gap is created such that a torque reversal should 
be accompanied by a minimum in the accretion luminosity. 
The model predicts a change in pulse profile from "leading- 
edge bright" during the spin down to "trailing-edge bright" 
during the spin up. Both, a minimum in luminosity and a 
chan ge in the shape of the pulse profile are observed in GX 
1 +4 dGreenhiU et al.l |l999), but neither effect is present in 
4U 1907-1-09 (see discussion above and Fig.^}. 

Recently, Perna et al. (2006) presented a new Ansatz to ex- 
plain torque changes in accreting neutron stars, which does not 
require the presence of retrograde disks. This model makes use 
of the fact that X-ray pulsars are oblique rotators, i.e., the neu- 
tron star's magnetic field is tilted by an angle x with respect 
to the axis of rotation of the neutron star. Assuming the accre- 
tion disk is situated in the neutron star's equatorial plane, for a 



ring of matter in the accretion disk, the magnetic field strength 
then depends on the azimuthal angle and thus th e boundary 
of the magnetosphere is asymmetric. As shown bv lPerna et alJ 
(2006), such a configuration can lead to regions in the disk 
where the propeller effect is locally at work, while accretion 
from other regions is not inhibited (i.e., on a ring on the disk, 
there are some regions with < r^o and other regions with 
'"ra > '"co)- This results in a nonlinear dependence of the accre- 
tion luminosity from the M through the outer parts of the disk. 
A nontrivial consequence of this more realistic accretion ge- 
ometry is that for values of x greater than a critical value, Xcnt, 
limit cycles are present, where cyclic torque reversal episodes 
are possible without a change in M. For typical parameters. 
;k'crit is between ~25° and ~45°. We stress that the existence of 
the limit cycles depends only on x ™d on the pulsar's polar 
magnetic field, B, and no variation of external parameters is re- 
quired to trigger torque reversals. This fact is a big advantage 
of this model over the models discussed above. 

Perna et al. (2006) show that their model can explain the 
principal observed properties of GX 1 +4, including the lumi- 
nosity drop during torque reversal, except for the observed cor- 
relation between torque and luminosity during the spin down 
phase. For 4U 1626-67, assuming B = 2.5 x lO'^G and;^' = 
68°, all observed properties of the spin history of 4U 1626-67 
including the large and values of \P/P\ before and after the 
torque reversal and the virtually unchanged luminosity of the 
source can be explained (in the model, the required change in 
luminosity is only ~5%). For these model parameters, the limit 
cycle of 4U 1626-67 is predicted to have a long period. 

The similarity between the pulse histories of 4U 1907 H-09 
and 4U 1626-67 suggests that a similar model would also work 
for the former source. We stress, however, that further quanti- 
tative work, which is outside the scope of this paper, is clearly 
required. Specifi cally, from a theo retical point of view we note 
that the model of IPerna et al ] 12006) does not yet take changes 
in M due to the orbital eccentricity into account and that the 
current version of the model assumes the disk to be flat and 
in the equatorial plane of the neutron star, while more realistic 
disks onto magnetized objects can be exp ected to be warped 
and precessing (e.g., Pfeiff^er & La3 l2004l and therein). From 
an observational point of view, the donor stars in both sys- 
tems are very different, with 4U 1626-67 being a low mass 
and 4U 1907-1-09 being a high mass system. Finally, we also 
emphasize that the torque reversal in 4U 1626-67 was from a 
spin up to a spin down, while for 4U 1907-1-09 it was from a 
spin down to a spin up. Since the torque reversal episodes are 
cyclic, however, we do not expect this latter difference to be of 
major importance. 

In conclusion, the torque reversal of 4U 1907 H-09 with no 
associated drop in luminosity and no change in the shape of 
the pulse profile seem to be difficult to reconcile with mod- 
els explaining the large P/P through the presen ce of a ret- 
rograde disk. On the other hand, the proposal of IPerna et alJ 
(.2006;) to explain torque reversals through accretion onto an 
oblique rotator seems to be able to explain the long term con- 
stancy of P trends in 4U 1907 -h09 and 4U 1626-67, the large 
magnitude of \P/P\, and the constancy of X-ray flux over the 
torque reversal, although further and more detailed theoretical 
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studies are required. As a prediction of the model for systems 
such as 4U 1907+09 and 4U 1626-67 is that torque reversal 
episodes are rare events, further monitoring of the pulse period 
of 4U 1907+09 is required to determine whether this change 
is a sign of a long term torque reversal or whether we are just 
observing a short deviation from a long term spin down trend. 
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